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ABSTRACT 


The  purpose  of  fhis  report  is  twofold:  (1)  To  review  some  aspects  of  the  theory  on 
density  fluctuations  in  a  plasma  (Port  I)  and  (2)  to  introduce  a  novel  method  for 
the  determination  of  the  fluctuation  spectrum  of  a  multicomponent  plasma  (Part  II). 
Part  I  attempts  to  follow  the  theory  that  starts  from  first  principles  (dynamics  of 
point-charged  particles)  and  leads  to  the  derivation  of  coarse-grained  properties 
of  a  plasma  in  order  to  compute  the  response  of  the  plasma  to  an  externally  ap¬ 
plied  longitudinal  field.  By  applying  the  Nyquist  theorem  to  the  response  func¬ 
tion  of  the  system,  we  derive  the  fluctuation  spectrum. 

Part  II  introduces  a  technique  by  which  the  spoce  and  time-dispersive  response 
function  of  a  plasma  that  has  an  arbitrary  number  of  charged  constituents  con  be 
evaluated  with  little  effort.  The  technique  consists  of  reducing  the  problem  to 
one  in  network  theory.  The  network  corresponding  to  a  plasma  is  a  bank  of  ca¬ 
pacitors  in  parallel,with  each  capacitor  representing  one  constituent  of  the  plasma. 
A  plier's  entry  at  any  desired  location  in  the  network  determines  the  response 
function  to  an  applied  external  potential. 
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INCOHERENT  SCATTERING  FROM  A  PLASMA* 


I.  THE  ELECTRONIC  SPECTRUM 


A.  INTRODUCTION 

1  -3 

Recent  work  on  incoherent  scattering  from  a  plasma  has  centered  on  the  calculation  of 
the  ionic  effect  on  both  the  spectral  quality  of  the  scattering  cross  section  and  the  integrated 
scattering  cross  section.  Since  the  cross  section  is  proportional  to  the  mean  square  of  the  fluc¬ 
tuation  of  the  density  of  the  system,  this  work  effectively  consisted  in  generalizing  earlier  re- 
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suits  obtained  for  the  fluctuations  in  an  electron  plasma  by  Pines  and  Bohm.  This  generalization 
demonstrated  that  the  contribution  to  the  cross  section  from  the  ion-carried  electron  fluctuations 
is  far  more  important  than  the  contribution  from  the  pure  electron  fluctuations.  However,  we 
believe  that  the  basic  issue  of  what  it  is  that  we  are  scattering  from  can  be  discussed  without 
carrying  along  the  cumbersome  formulas  that  include  the  effect  of  the  ions. 

In  this  report  we  attempt  to  follow  the  theory  that  starts  from  first  principles  (dynamics  of 
point  charged  particles)  and  leads  to  the  derivation  of  coarse-grained  properties  of  a  plasma 
(the  plasma  as  a  continuum)  in  order  to  compute  the  response  of  the  plasma  to  an  externally  ap¬ 
plied  longitudinal  field.  We  then  compute  the  energy  input  from  the  external  field  to  the  plasma 
and  show  that  it  is  nonzero,  even  in  a  plasma  without  collisions.  However,  when  there  is  a  dis¬ 
sipation  process  in  a  system,  a  fluctuation  is  associated  with  it.  By  applying  the  Nyquist  theorem 
to  the  system  we  proceed  to  derive  these  fluctuations.  We  further  indicate  that,  for  large  wave¬ 
lengths,  these  fluctuations  become,  in  effect,  organized  oscillations  at  the  plasma  frequency  a>p. 
We  then  conclude  that  in  a  phenomenological  way,  scattering  of  an  electromagnetic  wave  from  a 
plasma  can  be  described  as  scattering  from  a  density  wave. 

B.  THE  EQUATION  OF  MOTION  FOR  A  PLASMA 

A  complete  investigation  of  density  fluctuations  in  a  plasma  would  require  study  along  the 
5  6 

lines  described  by  Tchen  or  Rostoker  and  Rosenbluth.  Although  there  has  never  been  any  doubt 
as  to  the  soundness  of  the  basic  dynamic  principles  from  which  we  must  start  in  order  to  describe 
a  plasma  and  its  fluctuations,  the  statistical  processing  of  the  basic  equations  pertaining  to  many 
interacting  particles  (to  derive  macroscopic  properties  of  the  system)  introduces  considerable 
difficulties.  Many  compromises  and  mental  models  have  to  be  used  in  order  to  keep  a  clear  and 
useful  link  between  first  principles  (dynamics,  statistics)  and  final  results  (macroscopic  vari¬ 
ables).  By  keeping  track  of  the  meaning  of  the  compromises  and  of  the  models  we  may  hopefully 
retain  the  connection  between  a  microscopic  model  and  a  macroscopic  result.  To  this  end  we 
have  adopted  the  method  of  Klimontovich^  which,  though  deceivingly  simple,  is  very  helpful  for 
following  the  mathematical  manipulations. 

*  This  report  contains  notes  on  some  aspects  of  the  problem  of  incoherent  scattering  from  a  plasma.  These  notes 
were  prepared  in  the  summer  of  1961. 
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Consider  now  a  system  consisting  of  N  electrons.  When  electromagnetic  interactions  are 
neglected,  the  Hamiltonian  of  the  system  is 

N  N 

E  Pi^  +  U(q.))  +  |  E  V(q.-q.)  ,  (1) 

i=l 

—  —  th 

where  and  q^  are  the  momentum  and  the  position  of  the  i  "  particle,  U(qj)  is  an  external  po¬ 
tential  arising  from  an  external  force  acting  on  each  individual  particle  and  V{q,  —  q.)  = 

2  / 1—  —  I  ^ 

®  /  1 I  is  the  electrostatic  interaction  potential  of  the  particles  taken  in  pairs.  Introducing 
the  function  N(pqt), 

N 

N(pqt)  =  Z  *(p-Pj)  *(q- -  (2) 

i=l 

for  the  number  of  particles  at  pq,  we  can  write  for  any  function  g(pq), 

N 

j*  g(pq)  N(pqt)  dpdq  =  ^  g(Pj,qi)  . 

i=l 

Accordingly,  the  Hamiltonian  can  be  written  as 

^  ^  p^N(pqt)  dpdq  +  4  C  — ^ - ;  N(pqt)  N(p'q't)  dpdp'dqdq' 

+  U(q)  N(pqt)  dpdq  .  (3) 

We  now  note  that,  in  a  strictly  formal  way,  the  Hamiltonian  density  in  phase  space  can  be  written 
as 

AH  =  ^  p^N(pqt)  +  -4-  C  ;  N(pqt)  N(p'q't)  dp'dq'  +  U(q)  N(pqt)  ; 

|q-q-| 

therefore,  its  derivative  with  respect  to  N  is 

*=  =  =  4  p"  + 1  N(p'q't)  dp'dq-  +  U(q)  . 

Here  K  represents  the  Hamiltonian  of  one  particle;  therefore,  we  can  use  it  to  write  the  equa¬ 
tions  of  motion  for  this  particle. 


r  ~  —  N(p'q't)  dp'dq'  - 

3q  3q  |q  —  q'l  ®q 


where  p  and  q  are  the  canonical  variables.  Introducing  these  variables  in  the  continuity  equa¬ 
tion 
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3N(pqt)  ^  3N(pqt)  ^  ^  3N(pqt)  _  q 

aq  ap 

we  finally  obtain 

+  P  g^lpql)  -  (-1:  r  - .  N(p'q't)  dp'dq')  g-N(|q^l 

aq  'aq  |q  — q'l  '  9p 


_  auu)  aN(pqt)  ^  0  ^ 

aq  ap 

Equation  (4)  is  an  exact  equation  describing  the  microscopic  evolution  of  the  random  function 
— —  8 

N(pqt).  Gross  has  pointed  out  that  the  procedure  used  by  Klimentovich  poses  a  number  of  ques¬ 
tions;  nevertheless,  Gross  was  able  to  derive  the  same  equation  for  N(pqt)  by  using  a  more  di¬ 
rect  method.  In  any  case,  the  surprising  fact  is  that  we  can  obtain  an  exact  microscopic  equa¬ 
tion  for  the  evolution  of  N(pqt)  which  looks  so  much  like  the  well-known  transport  equation  of 
Boltzmann- Vlasov.  To  be  of  any  practical  use,  Eq.  (4)  has  to  be  averaged  in  order  to  derive  a 
macroscopically  meaningful  density  from  the  highly  discontinuous  N(pqt).  In  the  macroscopic 
description  of  the  system  we  are  not  interested  in  the  complete  set  of  mechanical  variables,  but 
only  in  a  much  more  restricted  number  of  variables  such  as  the  energy,  electric  currents  or 
charges  pertaining  to  macroscopically  infinitesimal  regions  of  the  system.  In  effect,  we  are  in¬ 
terested  in  a  formulation  in  which  density  of  matter,  charge  density,  etc.,  are  continuous  func¬ 
tions  of  position  and  time.  The  "homogenization"  of  the  discrete  system  is  made  through  the 
classical  statistical  methods  by  considering  an  ensemble  of  similar  systems.^  For  our  purpose 
we  need  only  note  that 

<N{pqt)>^^  =  f(pqt)  , 

where  f(pqt)  is  the  distribution  function,  usually  known  as  the  one-particle  momentum  distribu¬ 
tion  function.  With  this  step  the  individuality  of  the  particles  in  the  system  has  been  lost,  and 
we  must  begin  to  think  in  terms  of  a  fluid.  Consider  the  ensemble  average  of  all  terms  in  Eq.  (4), 

af^)  ^  _p  .  gflpjt)  _J_C  _e\  /N(p.q.t)  gN(pqtK 

aq  dq  ''  |q— q'l  8p  “V 


^  ^aU(q)  aN(pqtk  ^ 

'  aq  ap  '  av 

The  average  inside  the  integral  is  a  pair  distribution  function 
<N(pqt)  N(p'q't)>^^  =  f^Cpp',  qq',  t) 

denoting  the  extent  of  the  correlation  at  one  time  t  of  two  "particles, "  one  at  q,  the  other  at 
q'  and  with  respective  momenta  p  and  p'.  As  a  first  approximation,  we  assume  that  the  pair 
distribution  function  can  be  separated  as  (Ref.  7) 


f2(pp',  qq',t)  =  f(pqt)  f(p'q't)  +  ((i(pp',  qq',  t)  ,  (6) 

where  again  ^  is  a  pair  distribution  function.  It  can  be  argued  that  since  au(q)/dq  represents 
an  external  force,  not  necessarily  of  a  statistical  nature. 
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^  4  <N(5qt)>  =  -  •  4  f{Mt) 

'9q  8p  '  3q  3p  9p 

Introducing  Eqs.  (6'  and  (7)  into  Eq.  (5),  we  obtain 


(7) 


+  R  3£M  _  1.  r  _eL_  .  3f(pqt) 

9q  3q  I  q  —  q'  I  3p 


r  3;/»(pp'^qq'.t)  ,, 

9p  3q  I  q  —  q'  I  9p 

It  is  now  easy  to  obtain  the  standard  form  of  the  Boltzmann- Vlasov  equation  from  Eq.  (8),  and 

9 

also  the  equation  which  has  been  used  by  Gasiorowicz,  et  al.,  for  the  test  particle.  First  con¬ 
sider  F®^^(qt)  =  0  and  ^(pp'qq't)  -  0  and  call 

pintj-t)  =  _  4  C  f(p’q't)  dp'dq'  , 

9q  iq-q'l 


then 


9t  ^  m 


+  F 


9q 


int  ^ 

3p 


=  0 


(9) 


In  the  evaluation  of  the  integral,  q  =  q*  is  excluded;  therefore,  F^”^(  qt)  is  the  internal  force  act¬ 
ing  at  a  point  q  and  is  due  to  all  other  particles  in  the  system.  We  note  that  F^'^^  is  the  self- 
consistent  force  in  the  system  and  that  Eq.  (9)  is  an  approximate  version  of  Eq.  (5)  which  does 
not  take  account  of  the  pair  distribution  function.  We  know,  however,  that  the  two-particle  dis¬ 
tribution  function  (itself  a  function  of  the  one-  and  three-particle  distribution  function)  has  im- 
10 

portant  properties  leading  to  the  Debye  screening  potential.  It  is  therefore  expected  that  the 
simplification  leading  to  Eq.  (9)  would  require  further  justification.  Assuming  that  such  a  justi¬ 
fication  exists,  it  can  be  extended  to  the  case  where  ^  0.  Equation  (8)  then  becomes 


+  1.  MM)  _  A  r  dp'dq'  f(Pq.t)  MM) 

^  9q  Oq-J  Iq-q'l  3p 


=  -F^^(qt)  •  Mpqt)  ,  (10) 

3p 

This  is  the  equation  used  by  Gasiorowicz.  It  will  be  used  here  for  the  determination  of  the  re¬ 
sponse  of  the  system  to  an  applied  external  force  rather  than  for  the  response  of  the  system  to 
a  test  particle.  Since  the  observable  quantity  in  a  plasma  is  the  macroscopic  density  and  not  the 
distribution  function,  the  momentum  will  have  to  be  integrated  out  of  the  equations  to  determine 
the  density 

n( qt)  =  y  f(pqt)  dp  . 

Consider  then,  that  the  variations  of  the  density  from  the  average  value  are  due  to  variations  of 
the  distribution  function  from  its  equilibrium  value  fQ(p^), 

f(pqt)  =fQ(p^)  +f^(pqt)  . 
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(11) 


Introducing  the  above  into  Eq.  (10)  and  retaining  only  linear  terms,  we  obtain’' 


3f 


__1 

Sq 


8q 


“  r  !  ^  tTTix  V  ^ 


eext,-,., 


dp'dq'  .  f,(p'q't)  ^  =  _  F'’’^''(qt)  •  — 

|q-q'l  9p  8p 


Equations  (9)  and  (11)  are  widely  used  equations  for  the  distribution  function  f.  They  have 
deficiencies,  but  they  represent  an  improved  description  of  the  plasma  when  compared  with  the 
classical  Boltzmann  equation  with  the  collision  term.  The  main  reason  for  this  is  that,  in  an 
ionized  gas,  the  particles  are  constantly  under  the  electrostatic  field  set  up  by  all  the  othei  par¬ 
ticles;  in  a  neutral  gas,  the  particles  are  virtually  moving  in  free  space,  except  when  they  come 
close  to  each  other.  In  an  ionized  gas  "collisions"  are  not  discrete  events  unless  the  collisions 
are  with  neutral  particles  that  may  happen  to  be  in  the  plasma.  Whatever  success  has  been  de¬ 
rived  in  the  discussion  of  fully  ionized  gases  can  be  attributed  directly  to  the  electrostatic  term 
in  Eq.  (11).  Despite  its  shortcomings,  Eq.  (11)  can  at  least  be  used  to  derive  a  number  of  useful 
macroscopic  parameters  such  as  the  charge  density,  the  current  density,  and  the  electric  fields. 
In  the  following  we  shall  compute  the  charge  density  induced  in  the  system  by  an  external  poten¬ 
tial.  We  can  picture,  for  example,  a  plasma  confined  between  the  two  plates  of  a  capacitor  and 
the  external  field  due  to  a  battery  connected  to  the  plates.  For  all  practical  purposes  we  shall 
then  deal  with  a  dielectric  as  described  in  electrostatics. 


C.  THE  LINEAR  RESPONSE  TO  AN  APPLIED  FIELD 

The  density  induced  by  an  external  force  can  be  computed  in  a  straightforward  manner  if 
we  Fourier-transform  Eq.  (11).  Let 


^^(Pqt)  fl(p^^“)exp[ik  •  q-iwt]  — ~ 


(27r)" 


f^(pkw)  =  ^  f^(pqt)  exp[-ik  •  q  +  iwt]  dqdt 


Multiplying  Eq.  (11)  with  exp[-il{  •  q  +  iwt],  integrating  over  q  and  t,  and  noting  that  with 

q  -  q'  =  Q 

V(k)  =  r  exp  [-ik  •  Q)  V(Q)  dQ  =  r  exp  [-ik  •  Q]  dQ  =  , 

iQl  ^ 

we  obtain 

(  — iw  +  ik  •  p/m)  f .  (pkoj)  —  ik  \  f.  (p'q'w)  exp  [— ik  •  q']  dp'dq'-  — ^ 

k  ''  8p 

=  -F®’^(k«)  •  , 

8p 

or,  since  / f^(p'q'ai)  exp  [-ik  •  q*]  dp'dq’  =  5n(kw),  it  follows  that 


*The  term  inside  the  integral  with  Fq  has  been  neglected  since  it  is  canceled  by  the  neutralizing  positive  chorge 
in  the  system. 
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(  — ico  +  ik  • 

Dividing  by  (  — io)  +  ik  • 

6n(kw)  —  6n(kco)  ik 
or,  solving  for  6n(kco), 

£n( ku) 

.  ,  tjre  r  •  -  •  u 

T2  •'  — r  -r 

k  <0  —  k  •  p/m 

Equation  (12)  gives  the  desired  expression  for  the  response  of  the  system  to  external  force.  If 
this  force  is  of  electric  origin,  then  we  can  write  F  =  — eE  =  eV_^  or,  in  koi-space,  F(kai)  = 
eik^(kai).  The  induced  electric  charge  p(ku)  =  — e6n(ka))  can  then  be  written  as 

p(kcj)  =  C(ku)  <p®^(ka))  ,  (13) 

where  the  response  function  C(ka))  has  the  form 


2  9i  of 

p/m)  f . (pkoj)  —  ik  6n(kw)  — -  =  —  F®^^(kw)  •  -*-2 

^  k‘^  3p  ‘'P 

p/m)  and  integrating  the  resulting  expression  with  respect  to  p,  we  have 

2  ^  (3/3p)fo  _  n  (3/3p)fo 

"  dp  = -F®^^(ka))  ' 


4ire 


r  Ko/op)  Iq 
J  —  ico  +  ik  •  p. 


p/n 


r  wopMo 

—  io)  +  ik  •  p, 


p/n 


dp 


-ext-  r  (9/9p)fo  - 

,„ext,,...,  .  f  dp 


-iF^’^^ko;)  •  / 


fa)  —  k  •  p/n 


k  •  (3/9p)  f„ 


(12) 


C(kfa))  = 


k  •  (3/3p)  fp 
fa)  —  k  •  p/m 


dp 


.  .  4xe2  .k-  (9/8p)fo 
1  +  — y-  J  - -  dp 


fa)  —  k  •  p/m 


(14) 


The  use  of  C  has  been  introduced  in  order  to  indicate  that  what  we  have  computed  is  a  capacitance. 

This  point  underscores  the  fact  that,  at  this  stage  of  our  calculations,  we  are  dealing  with  a  model 

of  a  plasma  that  is  completely  divorced  from  any  concept  of  granularity  in  the  microscopic  scale. 

Before  discussing  the  properties  of  the  response  function,  we  must  proceed  with  computations 

that  will  give  a  more  appealing  form  to  C(kfa)).  For  this  purpose  we  adopt  the  notation  that  has 
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been  introduced  for  the  first  time  by  Lindhard.  This  notation  has  the  advantage  of  retaining 
some  of  the  concepts  of  electrostatics  in  continuous  media.  Lindhard  introduces  the  space  and 
time-dispersive  dielectric  constant  €(kfa))  in  the  following  way.  He  begins  with  the  linearized 
version  of  Eq.  (9) 


3t 


+  £■•—  +  F^"*(qt)  .  —  =  0 
3q  3p 


(15) 


which  he  proceeds  to  solve  along  the  lines  used  earlier  in  this  section. 
Eq.  (15)  and  solving  for  f^(pkfa)),  he  finds 


f^(pkfa))  =-iF‘'*^(kfa))  • 


0/3p)  fo 
fa)  —  k  •  p/m 


Fourier  -transforming 


therefore. 
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However, 


and 


_int  -  C  ^0  - 

6n(ka>)  = -iF''"^{ka>)  •  \  - r — z —  dp 

to  —  k  •  p/m 


(16a) 


F‘”*(kto)  =  eik(p^"^(kto) 


5p(kto)  =  — e6n(kco) 


P  k  ■  0/ap)  fo  _i  _ 

=  -e  \  - = - dp  <(7^"’^(kto) 

CO  —  k  •  p/m 


(16b) 


Let  us  now  turn  to  the  basic  theory  of  electrostatics  in  continuous  media.  We  know  from 
12 

this  theory  that: 

(a)  The  dielectric  polarization  in  the  medium  is  defined  by 

ind  ,  .  „ 

where  <p>  is  the  average  charge  in  the  medium. 

(b)  By  definition  the  internal  electric  field  is  given  by 

V  •  E  =  47r<p>  =  — 4jrV  •  P 

(c)  By  definition 

E  +  4jrP  =  D  . 

(d)  The  constitutive  equation  between  D  and  E  is  given  by 

D  =  eE  . 

From  the  above  it  follows  readily  that  in  kco -space  we  have 

p‘"^(k(o)  =  -ik  •  P(k(o)  =  -  ^  [D(kw)  -E(kto))  . 

On  the  other  hand  we  know  that,  in  a  frequency-dispersive  medium,  D(w)  =  €(to)  E(eo).  General¬ 
izing  this  relationship  to  a  space-dispersive  medium  we  write  D(]<co)  =  e(kto)  E(k(o),  therefore, 

p‘"'^(k(o)  =-^-  fe(kc^)-l]E{kw) 


=  -4^  [€(k&)) -1]  v’(kw) 


(17) 


Comparing  Eqs.  (17)  and  (16b),  we  then  have,  according  to  Lindhard,  the  definition  of  the  "longi¬ 
tudinal  dielectric  constant."  It  is  called  longitudinal  because  it  describes  electrostatic  effects  in¬ 
side  the  medium  and  is  given  by  the  expression 


€^(ka;)  =  1  +  ° 


JK  •  IV/ vp;  _ 

- -  dp 

o)  —  k  •  p/m 


(18) 


Returning  now  to  our  equation  for  the  induced  charge  due  to  an  external  potential  [Eq.  (13)] 
we  obtain 
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ext/T  V 
(p  ( ko)) 


(19a) 


and 


p(kaj)  = 


[€^(ka))  -  ll 
i  e^(ka>)  * 

I  e^(kco)J 


ext,r  \ 

<p  (kaj) 


(19b) 


C(ku)  =-^  [l --p4r-l  •  (20) 

^  I  e'"(kco)J 

For  the  explicit  determination  of  e^(ku))  we  need  only  note  that,  in  the  integral  of  Eq.  (18),  we 
have  to  make  the  following  interpretation  of  the  singularity; 

- ^ —  +  Y  ^(“  —  k  •  p/m)  , 

u  —  k  •  p/m  w  —  k  •  p/m 

where  ifi  denotes  the  principal  part.  Simple  integrations  then  give 


"(kcj)  =  1  + 


IR(ka))  +  il(ku))] 


where  is  the  Debye  length 


47re^n„ 


(21) 


and 


where  v^  =  ZwT/m  is  the  mean  thermal  velocity  of  the  electrons  at  temperature  T.  It  can  be 
easily  shown  that 


T  jW  X 

eL(kw)-.l-  (J) 

for 

J£_  »  1 

(24) 

for 

A 

A 

(25) 

It  is  important  that  at  the  limit  of  small  wave  numbers  the  longitudinal  dielectric  constant 
goes  to  tlie  classical  dielectric  constant  €(a))  since  we  find,  in  both  constants,  the  plasma  fre¬ 
quency  u)  .  However,  the  two  constants  have  a  completely  different  meaning  and  cannot  be  used 
P 

indiscriminately.  The  longitudinal  constant  implies  density  waves  similar  to  acoustic  waves 
with  a  frequency  Wp  and  a  wave  number  k  along  the  direction  of  propagation.  On  the  other  hand, 
the  classical  dielectric  constant  implies  a  propagation  of  a  transverse  electromagnetic  wave.  A 
space-dispersive,  transverse  dielectric  constant  can  be  developed  for  transverse  propagation  in 
order  to  generalize  the  classical  dielectric  constant,  but  it  is  outside  the  scope  of  this  report. 
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Let  us  use  now  the  two  limiting  forms  of  e^fku)  in  Eq.  (19).  For  a;/ kv^  »  1,  we  have 


P(k“)  -*  ^ 


2  2 
CO  —Cl) 

p 


^  D(w) 


(26) 


or,  using  the  continuity  equation. 


2  n 

CO  1 

_E _ 

2  2 

-“n  J 


E^’^^kco) 


For  co/kv^  «  1,  we  have 


p(kco)  -  - 


4Tr 


1 


1  +  k  X 


D 


ext.r 


272 


(27) 


The  first  limit  clearly  indicates  the  singularity  at  the  plasma  frequency,  the  second  limit  indi¬ 
cates  the  independence  of  the  response  function  from  the  frequency.  To  obtain  more  insight  we 
will  plot  the  function  (47r/k^)  C(kco)  =  D(co).  The  imaginary  part  of  D(co)  that  was  lost  while  the 

limit  for  co/kv  »  1  was  taken  can  be  reintroduced  by  the  substitution 
m 


Thus 


J -  1  ^  ’T 


,  +  —  6((o  ±  CO  ) 

COiCO  CO±CO  1  P 

p  p 


D(co)  =  2  2  =  -2^  [co-co  -  CO  +  CO  ) 

CO  —  CO  I  p  p  J 


=  [(P - ^ - —  <P  — 1 -  +  6(u)  -  CO  )  —  j-  6(co  +  CO  )] 

2  [co-cOp  CO+COp  i  p  1  pj 

The  functions  D'((o)  and  D"((o)  are  plotted  in  Fig.  1.  In  a  real  situation  the  singularities  will  be 
considerably  smoother.  This  is  indicated  by  the  dotted  curves.  Before  ending  this  section,  we 
shall  derive  some  more  useful  formulas. 

|3-Z5-4I»2| 


Fig.  1.  The  real  and  imaginary  parts  of  the  res^nse  function 
are  similar  to  the  dispersion  relations  for  an  oscillator. 
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The  use  of  Poisson's  equation  was  inherent  in  the  previous  treatment.  In  fact,  it  is  con¬ 
tained  in  the  equation 


F^'^^{qt)  =-~  r  ■_  f.(p'q't)  dp'dq' 

3q |q-q'l 


which,  in  ku -space,  takes  the  form 

J  fl(plcw) 


dp 


(28) 


or 

k<)))  =  —  4ire5n(  k<^))  =  47rp^'^^(kco)  , 

which  is  the  Fourier  transform  of 

V^.p^'^^qt)  =-4irp‘"^qt)  .  (29) 

If,  however,  external  charges  are  introduced  in  the  medium,  Poisson's  equation  would  have  to 
be  written  as 

V^<p(qt)  =  -  4irp^"*{qt)  -  4irp®’‘*(  qt) 
or 

k^ip(ka))  =  4irp^*'*(ka))  +  4)rp®^^(ka))  .  (30) 

11 

Introducing  Eq.  (17)  for  the  internal  charge,  we  then  have 

k^€^(kw)  ^(ku)  =  4jrp®*^(kto)  .  (31) 

Equation  (31)  gives  the  response  of  the  internal  potential  to  an  external  charge.  Equation  (19) 
can  then  be  considered  the  dual  of  Eq.  (31),  since  the  first  corresponds  to  a  circuit  driven  by  a 
voltage  source  and  the  second  corresponds  to  a  circuit  driven  by  a  current  source. 

Combining  Eqs.  (17)  and  (19)  we  obtain 

.p‘"*(kw)  =  .  (32) 

e^(kw) 

Further,  since  +  ^7®^,  we  have 

tot,r  \  c^(ka))  +  l  ext/T  . 

<p  (kw)  =  — -  (p  (ku)  .  (33) 

€^(ku7) 

To  complete  this  section  we  include  the  Kramers-Kronig  relationships.  These  relations 
result  from  the  principle  of  causality  and  should  be  applicable  to  any  well-behaved  response 
function.  Although  these  relations  have  to  be  reworked  for  a  space-dispersive  medium,  it  is  be¬ 
lieved  that,  for  low  thermal  velocities  of  the  electrons,  as  compared  with  the  velocity  of  light, 
no  serious  modifications  need  be  expected.  Separating  the  response  function  C(kco)  into  its  real 
and  imaginary  parts 

C(ku)  =  C'(kw) +iC"(ka))  , 
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the  Kramers-Kronig  relations  are 


and 


C"(kw)  =  --P  r  du}  . 

'  *  TT  J  a>  “  w 

The  following  relations  can  also  be  derived: 


(34) 


(35) 


(36) 


where  the  plus  implies  t  >  0  and  the  minus  implies  t  <  0.  Using  Eqs.  (34)  and  (20),  for  w  =  0, 
we  find  the  useful 


C'(kO)  =  J 


However,  from  the  explicit  expression  for  C,  we  also  find 

1 


-Re 


e^(kO)  -  1 


e'^kO) 


IT  J  O) 


Im 


€^(ka)) 


do) 


(37) 


O.  ENERGY  INPUT  INTO  THE  SYSTEM 

With  the  available  results  we  can  now  proceed  to  determine  the  energy  input  to  a  plasma  from 
an  external  force.  In  the  usual  theory  of  plasma,  this  energy  input  can  occur  only  because  of  the 
collisions  of  electrons  with  neutral  particles  in  a  gas  that  is  not  fully  ionized.  In  the  simplified 
form  of  this  theory  the  collision  frequency  is  the  important  parameter.  This,  as  we  have  noted 
in  Sec.  I-B,  implies  that  the  electrons  are  moving  in  the  system  on  straight  trajectories  until 
they  encounter  a  neutral  particle.  In  a  fully  ionized  gas  the  collisions  have  to  be  interpreted  on 
a  continuous  basis,  and  we  may  expect  that  the  continuous  interaction  of  a  test  particle  with  the 
field  of  all  the  others  may  lead  to  more  pronounced  "deflections"  than  a  single  short-range  colli¬ 
sion.  In  any  case,  these  continuous  collisions  should  lead  to  the  loss  mechanism  inside  the  sys¬ 
tem.  We  now  place  our  plasma  thermodynamically  in  an  insulated  container  and  electrically 

z 

between  the  plates  of  a  condenser.  Then  we  go  back  to  Eq.  (11)  which  we  multiply  by  (1/2)  mv  = 
(l/2m)  p^  and  integrate  with  respect  to  v  or  p.  We  expand  f  into  the  equilibrium  f^  and  f^,  but 
this  time  retain  the  nonlinear  terms.  We  note  that  terms  containing  (9/9p)  Iq  are  odd,  and 
therefore  yield  zero  on  integration.  What  remains  is 


The  first  term  gives  the  rate  of  change  of  the  kinetic  energy  density  K,  the  second  term 
gives  the  divergence  of  kinetic  energy  current  vector  ,  while  the  other  two  terms  give 
—  F  •  f  (p/m)fj(pqt)dp.  However,  F=— eEand,  since  by  definition  the  average  current  in  the 
system  is  -  (e/m)  J  pf^  (pqt)  dp  =  y(qt),  we  finally  obtain 


The  interaction  energy,  therefore,  is  given  by  the  right-hand  side  of  Eq.  (38), 


9<£  :=ext,— -rind,—., 

=  -E  (qt)  •  3  (qt) 


-ff 


E®^*(ka))  ■  exp[i(k  +  k')  q  -  i(a)  +  a>’)  t]^^^  da^ 

<27rr  (2irr 


Integrating  over  time  t  and  space  q,  we  then  obtain 


ff 


Tt 


E®^‘(kca)  •7i’'^(-k,-co)  ^ 

(27r)^ 


The  continuity  equation 

a 


gf  +  V  j=0 

becomes  p{ku)  =  k  •  j  (ka))/oj  in  kw-space.  Since  E(ka))  =  -ik(/)(kaj),  we  can  write 


-O)) 


dkdo) 

(2ir)'^ 


and,  remembering  Eq.  (19)  for  the  response  to  an  external  potential, 

^  dtdq  =  ia)<p'^^^( kaj)  C(— k,  —a))  — ca) 

=  Cria.k®='^ka,)|2c=Mka>) 

(2n) 


=  1  Cr  iw|<p®^\ka.)|^  [C>;'(ka>)-C(kw)]  . 

{Zirr 


or,  finally. 


JJ  ^  l'P®’‘*(ka.)|^a)Im[C(ka))] 

From  Eq.  (20)  it  follows  that 

k^ 

Im[C(ka))]  =  Im 


dkda) 


€^(ka)) 


(39) 


(40) 


The  results  we  have  just  obtained  are  very  important.  In  the  first  place,  we  have  shown 
that,  as  expected,  the  continuous  interaction  of  the  electrons  through  the  electrostatic  forces 
leads  to  dissipation.  Then  we  have  indicated  a  specific  method  for  evaluating  this  dissipation. 
Finally,  we  have  put  all  this  in  a  very  appealing  form  through  the  longitudinal  dielectric  constant. 

However,  the  most  significant  property  that  can  be  deduced  from  the  above  discussion  comes 
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from  a  consideration  of  the  theory  of  Callen  and  Welton.  According  to  this  theory,  whenever  a 
dissipation  process  exists  in  a  system,  a  fluctuation  is  associated  with  it.  This  theory  can  be 
considered  as  the  generalization  of  the  older  theory  of  Nyquist.  We  shall  not  reproduce  here  the 
derivations  of  Callen  and  Welton.  Nevertheless,  a  very  short  account  of  some  of  their  results 
may  be  helpful.  They  show,  for  example,  that  dissipation  involved  in  viscous  drag  on  a  moving 
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airborne  particle  leads  to  a  fluctuation  in  position  and  velocity  of  the  particle.  They  show  that 
the  acoustic  radiation  resistance  of  a  small  oscillating  sphere  leads  to  a  formula  for  the  pressure 
and  density  fluctuation  in  a  gas.  They  show  that  radiation  resistance  to  acceleration  of  an  elec¬ 
trical  charge  leads  to  a  fluctuation  of  the  electric  field  that  is  equivalent  to  black-body  radiation. 

In  the  next  section  we  shall  proceed  with  the  calculation  of  the  density  fluctuation  in  a  plasma. 


E,  DENSITY  FLUCTUATIONS 

14 

To  derive  the  fluctuation  spectrum,  Callen  and  Greene  argue  as  follows.  They  first  sup¬ 
pose  that  the  system  is  under  the  influence  of  an  external  force  from  t  =  — «  to  the  time  t  =  0 
when  the  force  is  removed.  The  system  is  then  left  on  its  own.  Let  us  compute  the  transient 
response  of  the  system,  using  the  results  obtained  in  Sec.I-C. 

From  Eq.  (1 9)  we  have 

P(qt)  =  r  C(k£o)  (p^^^ikco)  expfik  •  q  —  iajt]  (41) 

^  (2ir)^ 


If  the  externally  applied  potential  was  constant  to  t  =  0  and,  acting  at  the  position  of  the  "particle" 
at  q  =  0, 


then 


ext,—., 

<p  ( qt)  = 


0 


t.$  0 
t  >  0 


Introducing  this  into  Eq.  (37),  we  obtain 

dke^-^ 


p(qt)  =  -  ^ 


(2.^) 


TT  jC(kO)  + 


in  J 


C(ka))  e 


-iwt 


•dwl  <(>Q 


From  Eq.  (36)  and  for  t  >  0,  we  have 

C(k0)  =  Xp  r  C(k.)  e!^"* 
in  J  w 

and  therefore  Eq.  (41)  can  also  be  written  as 


(42) 


p(qt)  =-y  cosoitdoil  Vo  ■  (43) 

Equation  (43)  describes  the  evolution  (decay)  of  the  system  after  the  driving  force  was  removed. 

Callen  and  Greene  further  argue  that,  since  the  autocorrelation  of  the  density  in  time  as 
given  by 

R(t)  =  <p(t)  p(t  +  t)>  =y  dp'W(p')  p'  <T,  p'|p>  (44a) 

(where  W(p')  dp'  is  the  probability  of  finding  p  in  the  range  p'  <  p  <  p'  +  dp',  and  <t,  p'|p>  is 
the  expectation  value  of  p  at  a  given  time  if  it  is  known  that  p  had  the  value  p'  at  t  seconds 
earlier)  and  since  the  expectation  value  of  p  is,  in  fact,  given  by  Eq.  (43),  there  is  a  clear  con¬ 
nection  between  the  autocorrelation  and  the  response  function.  The  above  results  should  be 
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generalized  to  incljide  spatial  variations  that  are  applicable  to  our  problem.  Instead  of  Eq.  (44a) 
we  must  consider 

R(Q,t)  =  <p{q,t)  p(q  + q,  t  +  T)>  =  y  dp'W(p')  p'  <Q,t;p'|p>  ,  (44b) 

where  <Qt,  p'  |  p>  is  the  expectation  value  of  p  at  a  given  time  and  space  point  if  it  is  known  that 
p  had  the  value  p'  at  t  seconds  earlier  and  at  Q  centimeters  farther  up.  By  using  the  above 
arguments  we  then  find,  according  to  Callen  and  Greene,  that 

=  ,  ,4s, 

where  k  is  the  Boltzmann  constant  and  T  the  temperature.  Equation  (45)  is  the  result  we  wanted. 
Now  we  shall  seek  the  spectral  form  of  Eq.  (45).  We  first  note  that 


R(Q,  T)  =  lim  ^  p(  q,  t;  TV)  p(  q  +  Q,  t  +  t;  TV)  dqdt 
T,  V 


where,  as  usual, 


Rm  ^  JJ  p(  q,  t;  TV)  dqqt  =  Um  ^  J  ,j,  3 


Then  we  denote  by  S(kaj)  the  following  expression 


S(kw)  =  lim  p*(kw:  VT)  p(ka);  VT) 

T  V  ^ 


=  dQdr  exp(-ik  •  Q  +  iwr]  lim  ^  p(q,t;  VT)  p( q  +  Q,  t  +  t;  VT) 

T,  V 


dqdt 


Thus 


R(Q,  T)  =  JJ  S(ku))  expfik  •  Q  -  iwT] 


(2jr)" 


(46) 


Comparing  Eq.  (45)  with  Eq.  (46)  and  observing  that 
C'(kw)  =  C'(k, -w)  (even) 

C"(ku)  =-C''(k. -w)  (odd) 

we  find 


S'(k(o)=-(Z)^^y‘^>  , 


(47) 
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where  (2)  is  to  be  included  only  if  w  >  0  is  considered.  Equation  (47)  expresses  the  connection 
of  the  spectral  properties  of  the  fluctuations  in  a  system  in  terms  of  the  imaginary  part  of  the 
response  function.  More  explicitly  for  the  plasma 


S'(ka))  =  (2)  ^  Im 

We  can  now  use  the  Kramers-Kronig  relation  [Eq.  (34)]  to  integrate  Eq.  (47)  and  find  the  "in¬ 
tegrated  spectrum" 


€^(ka)) 


hi 

4i 


(48) 


y  S'(ka))  do)  j*  dco  =-/cTC'(kO) 


=  (cT 


e^(kO)-  1 
e^(kO) 


1  +  l/k\^ 


(49) 


_  2  k^/cT 
®  ''O  ~2“ 


4ire‘'njj  1  +  k^Xj^ 


k^X  ^ 

2 

■  ®  ”o  ,  .  .  2,  2 


1  +  k"X 


D 


(50) 


We  note  that  we  have  derived  the  fluctuation  spectrum  of  the  charge  density.  Since  p*(k)  p(k)  = 
e“8n’*'(k)  in(k), 

|8n(k)|^  =  no-  -f  ^  .  (51) 

By  taking  limits  for  wavelengths  (k  =  Zir/X)  very  much  larger  than  the  Debye  length  and  very  much 


smaller  than  the  Debye  length,  we  obtain 

[  "o 

kXD»l 

(52) 

|5n(k)|^  = 

1 

1  .  2 

[  nQ(kXp)  , 

kXjj  «  1 

(53) 

The  first  result,  giving  a  fluctuation  equal  to  the  mean  number  of  electrons  per  unit  volume,  is 
not  surprising  since  this  is,  in  fact,  the  fluctuation  of  any  system  in  equilibrium  composed  of 
noniateracting  particles.  What  is  intriguing,  but  hardly  unexpected,  is  that  the  fluctuation  for 
the  long  wavelength  limit  is  considerably  reduced.  The  intuitive  interpretation  of  this  result  is 
that  the  plasma  (in  a  phenomenological  way)  is  like  a  very  tight  membrane  which  can  only  execute 
oscillations.  To  understand  this  model  better,  we  have  to  look  at  the  spectral  quality  of  the  den¬ 
sity  fluctuations.  The  spectrum  is  given  by  Eq.  (48)  (divided  by  e^).  We  have,  more  explicitly. 


|8n(ka))|  Vng  = 

4ire  n 


Im{e^(ku))) 
«■&)  |c  (ka))|^ 


(54) 


■  kX 


Fig.  2.  Mean-square  fluctuation  of  electron  density  in  a  positive  ion  sea. 


Fig.  3.  Gabor's  conception  of  the  dispersion  relation. 


since -Iin{e^' (kw)}  =  Im  {e^(ka))}.  With  the  help  of  Eqs.  (2a)  and  (23)  we  find  that,  in  the  limit 
of  w/kv^  »  1,  Eq.  (54)  becomes  sharply  peaked,  indicating  that  the  fluctuation  spectrum  loses 
the  quality  of  a  random  phenomenon.  In  fact,  it  does  not  represent  a  fluctuation,  but  an  oscilla¬ 
tion  at  a  frequency  very  close  to  the  plasma  frequency  co^. 

Figure  2  plots  the  mean-square  fluctuation  of  the  electron  density  (in  the  absence  of  ions). 

We  note  that,  as  an  observation  limits  itself  to  smaller  and  smaller  k's  or  larger  wavelengths, 
the  spectrum  of  |6n|  becomes  concentrated  more  and  more  near  one  frequency.  In  the  other 
limit  (of  small  wavelengths),  the  spectrum  is  Gaussian,  as  it  should  be  for  a  neutral  gas.  Smce 
the  sharp  peak  near  w  i  corresponds  to  a  longitudinal  oscillation,  we  have  obtained  the  natural 
motion  of  the  plasma  in  equilibrium. 

To  complete  this  picture,  we  shall  now  refer  to  the  analysis  of  Gabor^^  for  the  density  fluc¬ 
tuations.  First  we  note  that  the  dispersion  relation  u  =  a)(k)  for  the  plasma  can  be  obtained  from 
the  numerical  solution  of 


e^(kaj)  =  0 

which  gives  the  well-known  relation 


(55) 


co^(k)  i  k^ 

p  m 


(56) 


The  exact  solution  of  Eq.  (55)  would  give  a  curve  as  shown  by  the  solid  line  in  Fig.  3.^  Now  this 
curve  can  be  considered  as  the  projection  of  the  relative  maxima  of  Fig.  2  on  the  plane 
(cu/kvQ,  kXj^).  Referring  to  the  analysis  of  Gabor,  who  found  that  at  large  k's  the  dispersion 
curve  effectively  becomes  a  diffuse  line  describing  the  degeneration  of  organized  oscillations  to 
a  random  motion  of  the  particles,  we  may  plot  a  similar  curve  but  with  "equipotentials"  around 
the  classical  dispersion  relation  (Fig.  4).  Thus  the  dispersion  relation  is  a  diffuse  line  for  all 
k's  but  becomes  sharper  as  k  is  decreased.  Gabor  succinctly  expressed  this  behavior  of  the 
plasma  as  follows; 


"For  waves  much  shorter  than  the  Debye  length  there  will  be  no  dispersion  law, 
but  merely  a  [statistical)  correlation  between  frequency  and  wave  number." 

It  is  the  same  picture  that  we  get  from  the  work  of  Pines  and  Bohm^  who  find  that  the  density 
fluctuations  can  only  have  wave  numbers  smaller  than  kj^.  For  k  >  k^  the  plasma  behaves  sta¬ 
tistically  like  a  set  of  neutral  particles  interacting  with  short-range  forces, 

From  the  above  discussion  we  conclude  that  a  very  useful  model  of  scattering  from  a  plasma 
can  be  constructed  if  the  plasma  is  viewed  as  a  continuum  whose  density  varies  periodically  in 
space  and  time.  We  have  to  keep  in  mind,  however,  that  these  macroscopic  density  variations, 
or  fluctuations,  reflect  in  our  measuring  scale  the  granular  nature  of  matter  in  the  microscopic 
scale.  In  general,  this  granular  structure  is  not  accessible.  What  we  see  and  measure  are  its 
effects  on  our  measuring  apparatus.  In  the  particular  case  of  a  plasma,  we  can  conceive  of 
three  distinct  measuring  systems.  We  can  place  the  plasma  between  the  plates  of  a  capacitor 
and  measure  the  fluctuating  voltage  or  current.  We  can  send  an  electromagnetic  wave  through 
the  plasma  and  measure  scattering  due  to  fluctuations  in  the  density  (a  continuum  with  uniform 


t  In  terms  of  x  =  (mAvJ  and  a  =  kXp,  Eq,(56)  becomes  =  (2x^-3)"' . 
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density  does  not  scatter  radiation).  We  can,  finally,  send  a  beam  of  electrons  through  the  system 
and  measure  the  interaction  of  the  beam  with  the  plasma. 


Fig.  4.  Dispersion  relation  in  a  plasma. 
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II.  THE  ION- CARRIED  ELECTRONIC  SPECTRUM 


A.  INTRODUCTION 

In  Part  I,  we  gave  a  general  picture  of  the  models,  concepts,  and  techniques  used  to  establish 
the  fluctuation  spectrum  of  a  plasma  composed  of  electrons  in  a  neutralizing  sea  of  positive  charge. 
We  reached  the  conclusion  that,  on  a  macroscopic  scale,  the  fluctuation  spectrum  of  the  electrons 
could  be  interpreted  as  a  density  wave  propagating  in  the  plasma  with  a  given  frequency  to  and  a 
wave  number  k,  provided  the  measuring  scale  is  larger  than  the  Debye  length.  It  would  be  rea¬ 
sonable  to  expect  that  similar  results  could  be  obtained  for  a  model  plasma  in  which  the  role  of 
ions  and  electrons  is  reversed  (where  the  ions  are  the  particles  moving  in  a  neutralizing  negative 
charge).  The  only  difference  between  the  two  models  would  come  from  the  fact  that  ions  are 
heavier.  Since  the  characteristic  parameter  =  co/kv^  ^  in  the  fluctuation  spectrum  of  the  elec¬ 
trons  is  inversely  proportional  to  the  mean  thermal  velocity  v^  ^  =  JZKT/m^,  where  m^  is  the 
mass  of  the  electrons,  it  is  clear  that  the  results  obtained  in  Part  I  could  be  scaled  and  used  for 
the  ions  in  a  negative  continuum  if  a  new  parameter  were  to  be  defined  as 

X.  =  x_  \lM/m  =  {a)/kv„  )  VM/m  , 

where  M  is  the  mass  of  the  ions. 

In  reality,  however,  a  plasma  can  be  composed  of  electrons,  negative  ions,  positive  ions, 
and  neutrals:  therefore,  the  models  with  positive  or  negative  neutralizing  seas  can  only  be  use¬ 
ful  under  restricted  circumstances.  In  general,  it  is  expected  that  the  fluctuation  spectrum  of 
any  of  the  constituents  of  the  plasma  will  be  seriously  affected  by  the  presence  of  the  other  con¬ 
stituents.  On  an  intuitive  basis,  which  rests  on  the  assumption  of  perfect  neutrality  in  a  plasma, 
we  should  expect  that  ion  fluctuations  (or  ion  density  waves)  will  be  closely  followed  by  electron 
density  waves  or  fluctuations.  On  these  grounds  alone  and  because  the  ionic  mass  M  is  much 
larger  than  the  electronic  mass  m,  giving  an  ion  plasma  frequency  a>  .•  =  w  sIm/M,  we  can 
make  a  guess  that  the  portion  of  the  electronic  spectrum  between  w  =  0  and  w  =  Wp  will  definitely 
contain  an  important  contribution  due  to  ionic  waves.  The  problem  is  to  compute  the  exact  amount 
of  this  contribution.  As  stated  in  Part  I,  all  the  effort  in  the  theoretical  description  of  the  inco¬ 
herent  scattering  from  the  ionospheric  plasma  has  centered  in  the  evaluation  of  the  effect  of  the 
ions  on  the  spectrum  of  the  electron  fluctuations.  The  purpose  of  Part  11  is  to  introduce  a 
simple  technique  by  which  the  response  function  of  a  plasma  with  an  arbitrary  number  of  charged 
constituents  can  be  evaluated  with  little  effort.  Nevertheless,  we  shall  also  indicate  the  method 
by  which  the  equation  of  motion  for  electrons  [Eq.  (4)  of  Part  I]  can  be  modified  to  include  the  ef¬ 
fect  of  the  ions. 

B.  THE  EQUATION  OF  MOTION  FOR  THE  ION-ELECTRON  PLASMA 

The  generalization  of  Eq.  (4)  of  Part  I  to  derive  the  equation  of  motion  for  the  ion-electron 
plasma  can  be  obtained  by  considering  the  Hamiltonian  of  the  system  with  electrons  and  ions. 

This  Hamiltonian  contains  pair  Coulomb  interactions  between  electrons  and  electrons,  ions  and 
ions,  ions  and  electrons.  If  N  (pqt)  is  the  density  of  electrons  and  N^(pqt)  the  density  of  ions, 
from  the  procedure  described  in  Sec.  I-B,  we  can  easily  obtain  the  following  exact  pair  of  coupled 
equations. 
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aN  (pqt)  -  aN  (pqt)  ^  aN  (pit) 

- :rr -  +  -C.  -  —  —  F  Int)  -  +  P» 


-  —  E'=^‘'(qt) 
m 


aq  J  |q-q' 


X  {N.(p'q't)  -  Ng(p'q't)} 


2N^(pqt) 


aN.(Mt)  -  aN.(iit)  /  9N.(pit)\  2  a  f 

r  -  -  -  -  1  ®N.(pqt) 

x{N.(p'q't)-N  (p'q't)}  •  -  -  =0  . 

'  ®  9p 


where  E®  is  the  external  field.  The  parenthesis  in  Eq.  (58)  indicates  that  the  external  field 
acting  on  the  ions  can  be  disregarded  if  we  are  interested  in  the  response  of  the  electrons.  Since 
an  externally  applied  field  acts  on  both  electrons  and  ions,  it  is  obvious  that  disregarding  the 
term  in  parenthesis  implies  that  the  external  field  acting  on  electrons  readjusts  itself  to  a  differ¬ 
ent  value. 

The  pair  of  equations  for  and  N.  can  now  be  homogenized  in  order  to  obtain  the  coarse¬ 
grained  properties  of  the  ion-electron  plasma.  Writing 

<Ng(pqt)>^^  =  fg(Mt) 


<N.(pqt)>^^  =  f.{pqt) 


<N  (p'q't)  N  (pqt)>  =  <N  (p'q't)>  <N  (pqtj>  +  (p'q'pqt) 


^  fg  (p'q't)  fg(pqt)  + 


<N.(p'q't)  Ng(pqt)>g^  =  <N.(p'q't)>  <Ng(pqt)>  +  <;fjg(p'q'pqt) 


=  f. (p'q't)  fg(pqt)  + 


and  corresponding  equations  for 


<N.(p'q't)  N.(pqt)>g^^  ,  <Ng(p'q't)  N.(pqt)>^^ 

we  obtain  a  set  of  coupled  equations  by  expanding  the  single  particle  distribution  functions  in  the 
usual  manner  (f  .  =  f„  .(p)  +  f.  _  ;(pqt))  and  neglecting  all  the  pair  correlations  (p  ,  (p..,  (p  w- 

cc  IX  6X  Xli 


^  =  -5.  E®^(5t)  ^  (59) 

3p  ap 


3q  9q  |q-q'l 
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Following  the  same  technique  as  in  Sec..I-B,  we  can  determine 

and  also  any  desired  response  function.  From  this  point  on,  the  problem  becomes  a  problem 
of  algebra.  It  is  obvious  that,  if  more  components  are  introduced  in  the  plasma,  the  equations 
increase  in  number.  The  calculations  required  to  compute  the  induced  charges  or  the  response 
functions,  also  increase.  To  avoid  these  complexities,  we  have  developed  a  method  which  re¬ 
duces  the  problem  to  one  in  network  theory.  This  method  is  discussed  in  the  next  section. 


C.  RESPONSE  FUNCTION  OF  A  COMPOSITE  PLASMA 


In  Sec.  I-C  we  argued  that  the  plasma  can  be  described  in  terms  of  the  classical  theory  of 
dielectrics.  We  shall  now  show  that  this  model  of  a  plasma  can  lead  effortlessly  to  the  calcula¬ 
tion  of  the  response  function  of  any  constituent  of  the  plasma  in  the  presence  of  the  remaining 
constituents.  Consider  first  a  dielectric  on  which  we  have  placed  an  external  charge  The 

presence  of  the  external  charge  polarizes  the  dielectric,  or  in  other  words,  induces  a  charge 
inside  the  dielectric.  The  potential  inside  the  dielectric  is  given  by  Poisson's  equation 

V^<p(qt)  =  -  4irp^"‘*(qt)  -  4)rp®’^(qt)  (6la) 

or  in  ko) -space 

k^<p(kaj)  =  47rp^*'^(ka))  +  4)rp®^(ka))  .  (6lb) 

From  Eq.  (17)  of  Part  I, 


ind 

P 


[e^(k<o)  —1]  <D(kw) 


(62) 


relating  the  induced  charge  to  the  induced  potential  and  the  Poisson  equation,  we  can  then  obtain 

k^«^(ka))  ^)(ka))  =  47rp®*^(k&>)  ,  (63) 

the  generalization  of  Poisson's  equation  in  matter  and  for  a  space-dispersive  medium.  Suppose 
now  that  the  medium  consists  of  a  number  of  different  charged  particles.  It  is  evident  that  the 
external  charge  brought  to  the  dielectric  from  outside  will  induce  different  charges  on  each  of 
the  constituents.  If  each  of  the  constituents  were  to  be  described  by  a  dielectric  constant  of  its 
own,  then  we  could  write 

p.^"'^(kw)  =  -  ^  f€j^(k(o)  -1  ]  v’j(kw)  ,  (64) 

where  the  subscript  j  denotes  the  constituent.  It  is  clear,  however,  that  the  potential  in  the 
system  has  to  adjust  itself  so  that  it  will  be  common  to  all  constituents;  therefore,  <p.(kco)  =  <p(ku). 
Since  the  total  induced  charge  has  to  be  the  sum  of  the  individual  induced  charges,  we  can  write, 
for  the  Poisson  equation  in  kw -space. 


k^^(kai)  =  — k^  fc-^(kci»)  — IJ  ^(ku)  +  4)rp®*^(kw) 

i=i  ^ 


(65) 
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In  an  electrostatic  notation  we  could  further  write 


w'here 


C.(kco) 


) 


Co(k-)=-iF 


(66a) 


C.(kcj) 


[€^(ka))-l] 


j  ¥=  0 


(66b) 


The  above  discussion  suggests  that  the  composite  plasma  can  be  described  as  a  bank  of  "capaci¬ 
tors"  driven  by  a  charge  source.  The  equivalent  network  of  the  composite  plasma  is  shown  in 
Fig.  5.  This  corresponds  to  the  physical  situation  where,  for  example,  an  electron  is  shot  at  a 
volume  containing  a  plasma.  It  can  also  correspond  to  a  situation  in  which  neutral  molecules  in 
the  plasma  are  photo-ionized;  in  this  case  corresponds  to  one  electron  and  one  positive  ion 
beginning  their  existence  at  some  point  in  the  plasma  and  moving  with  relatively  high  velocities. 


Fig.  5.  The  equivalent  circuit  of  a  multicomponent 
plasma  driven  by  a  charge  source. 


The  above  description,  however,  is  not  adequate  for  the  study  of  fluctuations  in  any  of  the 
components  of  the  plasma.  We  have  seen  in  Sec.  I-C  that  for  this  problem  we  need  to  know  the 
induced  charge  on  the  component  in  question  when  an  external  potential  is  acting  on  this  compo¬ 
nent.  Had  the  network  of  Fig.  5  been  a  real  electric  network,  then  we  could  compute  the  desired 
response  by  making  a  "plier's  entry"  at  the  branch.  Introducing  a  potential  at  this  entry,  we 
could  calculate  the  charge  induced  on  the  capacitor,  and  therefore  the  response  function. 

This  is  equivalent  to  the  calculations  that  would  have  been  required  to  solve  the  system  of  equa¬ 
tions  of  Sec.  II-B.  The  network  for  the  calculation  of  the  response  to  an  applied  potential  is  shown 
in  Fig.  6.  It  is  obvious  that  we  can  use  this  network  for  the  calculation  of  any  desired  transfer 
function,  that  is,  the  response  at  the  i^^  capacitor  which  is  the  consequence  of  an  applied  potential 
through  a  plier's  entry  at  the  branch. 

I3-25-4H5  lfl-b)| 
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Fig.  6.  The  equivalent  circuit  of  a  multicomponent  Fig.  7.  The  equivalent  circuit  for:  (a)  electrons 
plasma  driven  by  an  external  potential.  in  a  positive  sea  and  (b)  electrons  and  ions. 

As  a  first  application,  consider  a  plasma  with  electrons  in  a  positive  sea.  The  equivalent 
network  is  shown  in  Fig.  7(a-b).  Since  the  capacitors  are  in  series,  we  have 
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^ext(k^)  _  ^®^^(kco) 


Introducing  the  values  of  the  C's  from  Eq.  (66b),  we  have 


Cg{kco)  =- 


j^2  e^hko,)-! 

4t!  L,r-  > 
Eg  (koj) 


(67) 


which  is  the  same  as  the  value  obtained  in  Part  I,  Eq.  (20).  The  second  application  concerns  the 
electron-ion  plasma.  From  the  circuit  of  Fig.  7(b)  we  obtain 


Cee(ka,)  = 


(Cq  +  C.) 

Cq+Cj  +  C^ 


(68) 


^2  e^iku)  [€^(ka))-l] 
€.^(kw)+  £^(ka))  —  1 


This  is  the  response  function  required  to  compute  the  fluctuation  spectrum  of  the  electrons  in  the 
presence  of  ions  and  is  pertinent  to  the  incoherent  scattering  from  the  ionospheric  plasma.  Pro¬ 
ceeding  with  the  same  technique,  we  find  that  the  response  of  a  plasma  with  different  types  of 
ions  is 


Cee(kco) 


kl 

4ir 


1  -r  +  (  S 
e^(ka>)  + 


ej^(ka)))  •  (e^ikco)  -  ij 


(jfl 


(70) 


where  r  is  the  number  of  types  of  ions.  In  a  similar  way,  we  can  determine  the  transfer  re¬ 
sponse  function  for  a  two-component  plasma  when  the  potential  is  applied  at  the  Cg  branch  and 
the  response  is  measured  at  the  C.  branch.  The  result  is 


<i<k“)=-fe 


[€^^(ka>)-l)[6.^(kw)-l] 
ej^(ka))  +  €^(kw)  -  1 


(71) 


We  can  also  find,  by  inspection,  the  total  response,  that  is,  the  total  charge  induced  when  the 
applied  potential  acts  on  and  Cj.  The  result  is 


k^  ^ 

Cg+i(kco)  =  ^  - 


^{ku)  +  e.^(ku))  —  2 

^(ko;)  +  €/"(k«)  -  1 


(72) 


This  completes  the  calculations  for  the  derivation  of  a  few  of  the  response  functions  required  to 
describe  a  plasma.  These  functions  can  be  put  to  immediate  use  for  the  calculation  of  the  mean- 
square  fluctuations  in  a  plasma.  Referring  to  Sec.  I-E,  we  find  that  the  spectral  quality  of  the 
fluctuations  is  given  by 

S(kw)  ^  Im  {-C(ka))} 


(73) 


where  C(kt<j)  can  be  any  of  the  functions  we  have  derived  in  this  section.  Since,  as  a  result  of 
the  Kramer-Kronig  relations,  the  integrated  spectrum  is  given  by  the  real  part  of  C(ka))  eval¬ 
uated  at  o)  =  0, 

y  S{ku)  da>  =  kT  Re  {-C(kO)}  ,  (74) 

we  can  compute  the  integrated  density  fluctuation  of  any  component  with  elementary  algebra.  In 
fact,  referring  again  to  Part  I,  Eq.  (25),  we  have 

€  .^(kO)  =  1  +  (1/kXj^.)^  ,  (75) 

where  now 


4rre^n. 

1 


In  a  multicomponent  plasma  we  must  have 


(76) 


Z  n.  =  n 
Li  j  o 

j=l 

so  that  the  total  number  of  the  different  kinds  of  ions  is  equal  to  the  number  of  electrons.  If 
negative  ions  are  to  be  included,  we  must  impose  the  condition. 


r  s 

E  nj(+)  =  %  +  Z  «j(-)  • 

3=1  3=1 

Using  the  above  results  in  Eq.  (69)  we  readily  obtain 


[«nee(k)( 

Te=T. 


1  +  k^Xp 

2  k\p 


(77) 


which  is  the  familiar  result  of  the  density  fluctuations  in  an  electron-ion  plasma.  With  equal  ease 
we  find  the  fluctuation  correlation  between  electrons  and  ions  from  the  transfer  response  function. 
For  equal  temperatures  it  has  the  value 


I  I 


°  2  +  2 


(78) 


Finally,  the  total  density  fluctuation  is 


dn 


tot ' 


2  +  k^X^ 


(79) 


Many  more  results  can  be  computed  (for  example,  |6n.(k)l^)  but  these  are  left  as  exercises  for 
the  reader.  The  spectral  shape  of  the  fluctuations  can  be  computed  from  the  response  functions, 
but  the  calculations  are  involved  because  of  the  many  variables.  In  Fig.  8  we  present  a  three- 
dimensional  plot  of  the  fluctuation  spectrum  of  the  ion-electron  plasma  when  the  temperatures 
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of  the  two  constituents  are  the  same.  This  plot  represents  the  value  of  S(ka))  in  Eq.  (73)  when 
C(kai)  is  given  by  Eq.  (69). 


Fig.  8.  The  fluctuation  spectrum  of  on  ion-electron  plasma  as  a  function 
of  X  =  o/kv0.  The  parameter  is  lAX))/  where  k  is  the  wave  number  and 
Xq  is  the  Debye  length. 
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